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We report a Dalitz-plot analysis of the charmless hadronic decays of charged B mesons to the final 
state K if^it . Using a sample of 226.0 ± 2.5 million BB pairs collected by the BABAR detector, we 
measure the magnitudes and phases of the intermediate resonant and nonresonant amplitudes for 
both charge conjugate decays. We present measurements of the corresponding branching fractions 
and their charge asymmetries that supersede those of previous BABAR analyses. We find the charge 
asymmetries to be consistent with zero. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The properties of the weak interaction, the complex 
quark couplings described in the Cabibbo-Kobayashi- 
Maskawa matrix (CKM) elements p| as well as models of 
hadronic decays, can all be studied through the decay of 
B mesons to a three body charmless final state. Studies 
of these decays can also help to clarify the nature of the 
resonances involved, not all of which are well understood. 
The decays B ± -> K ± tt^tt^ can proceed via intermedi- 
ate quasi two-body resonances as well as nonresonant de- 
cays. The interference among these resonant and nonres- 
onant decay modes can provide information on the weak 
(CP odd) and strong (CP even) phases. Theoretical pre- 
dictions using various hadronic decay models exist for the 
branching fractions and CP asymmetries of the decays 



nmetries or tne decays 

(770)^+ Huang. 



B+ -> K*°(892)7r+ and B+ -> p ( 
Precise measurements of the branching fractions of these 
modes and the CP asymmetry of B + — > p a (770)K + can 
discriminate among these models. The CP asymmetry 
of B + — > K*°(892)tt + is predicted to be zero, or at least 
very small, by all theoretical models. A measurement 
of a large asymmetry in this mode would therefore be a 
possible indication of new physics. 

In this paper we present results from a full ampli- 
tude analysis for B ± — > K ± tt t tt ± decay modes based on 
a 205.4fb _1 data sample containing 226.0 ±2.5 million 
BB pairs (N Bg -). These data were collected with the 
BABAR detector at the SLAC PEP-II asymmetric- 
energy e+e~ storage ring Q operating at the T(4S) res- 
onance with center-of-mass energy of y/s = 10.58 GeV. 
An additional total integrated luminosity of 16.1 fb _1 
was recorded 40 MeV below the T(4S) resonance and was 
used to study backgrounds from continuum production. 

A number of intermediate states contribute to the de- 
cay B^ — > K ± ir^TT :il . Their individual contributions are 
obtained from a maximum likelihood fit of the distribu- 
tion of events in the Dalitz plot formed from the two 
variables x — ro]L-± =f and y 



X 7T±7r=F ' 



Neglecting, at 



this stage, possible variations of the experimental accep- 



ts This version of the paper corrects a coding error discovered after 
submission to Phys. Rev. D and described in an Erratum submit- 
ted 26 Oct 2006. 



tance over the Dalitz plot, the probability density func- 
tion (PDF) for signal events is given, in the isobar for- 
malism (see for example [fj, [l(J, [ll| ) , by: 



P(x,y) 



I ILjCj-e^foy)! 2 dxdy' 



(1) 



The amplitude for a given decay mode j is Cjt j Fj(x, y) 
with magnitude Cj and phase 0j (—tt < 0j < tt). The 
magnitudes and phases are measured relative to one of 
the contributing channels, _ftT*°(892) in this analysis (c = 
1, 9 = 0). The distributions Fj describe the dynamics of 
the decay amplitudes and are a product of the invariant 
mass and angular functions. Examining the case where 
the resonance is formed in the x variable we have: 



Fj(x,y) = Rj(x) x Tj(x,y). 



The Fj are normalized such that: 



\Fj(x,y)\ dxdy = 1. 



(2) 



(3) 



For most resonances in this analysis the Rj are taken 
to be relativistic Breit- Wig ner lineshapes with Blatt- 
Weisskopf barrier factors |l2( • There is no cut-off applied 
to the Rj and so they are integrated over the entire Dalitz 
plot. 

The Breit-Wigner lincshape has the form 



Rj (m) = 



1 



(mg — to 2 ) — imoT(m) ' 



(4) 



where too is the nominal mass of the resonance, m is the 
mass at which the resonance is measured and r(m) is the 
mass-dependent width. In the general case of a spin J 
resonance, the latter can be expressed as 



r(m) = r 



2,7+1 



/TOo\ Xj(q) 
V m J X 2 j(qo) ' 



(5) 



The symbol Tq denotes the nominal width of the reso- 
nance. The values of too and Tq are obtained from stan- 
dard tables The value q is the momentum of either 
daughter in the rest frame of the resonance. The symbol 
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qo denotes the value of q when m = mo. X.j(q) represents 
the Blatt-Weisskopf barrier form factor |l2j : 

Xj =0 (z) = 1, (6) 
X J=l (z) = Vl/(l + z 2 ), (7) 
X J=2 (z) = ^17(^+3^+9), (8) 

where z = rq and r, the meson radius parameter is taken 
to be 4.0 GeV -1 Q. 

For the / (980) the Flatte form [TH| is used. In this 
case the mass-dependent width is given by the sum of 
the widths in the tttt and KK systems: 

r(m) = T v (m)+T K {m), (9) 

where 

IV(m) = gn\/m 2 - 4to 2 ( 10 ) 

T K (m) = g K \Jm 2 - Am 2 K (11) 

and and gx are the effective coupling constants 
squared for /o(980) -> tt+tt" and / (980) -> if+if - , re- 
spectively. Below the K + K~ threshold the function con- 
tinues analytically, the Yk term contributing to the real 
part of the denominator. 

For the angular distribution terms Tj we follow the 
Zemach tensor formalism 0, ^J. For the case of the 
decay of a spin B meson into a s pin J resonance and a 
spin bachelor particle this gives |18| : 

rpj=0 _ -, 
3 ~ 

T/ =1 = -2p-q (12) 

where p is the momentum of the bachelor particle and 
q is the momentum of the resonance daughter with the 
same charge as the bachelor particle, both measured in 
the rest frame of the resonance. 

The B candidates are reconstructed from events that 
have four or more charged tracks. Each track is required 
to be well measured and to originate from the beam spot. 
The B candidates are formed from thrce-charged-track 
combinations and particle identification criteria are ap- 
plied to reject leptons and to separate kaons and pions. 
The average selection efficiency for kaons in our final 
state that have passed the tracking requirements is about 
80% including geometrical acceptance, while the average 
misidentification probability of pions as kaons is about 
2%. 

Two kinematic variables are used to identify signal B 
decays. The first variable is AE = E* B — yfs/2, the differ- 
ence between the reconstructed center-of-mass (CM) en- 
ergy of the B- meson candidate and \fs/2, where yfs is the 
total CM energy. The second is the energy-substituted 
mass m ES = yj (s/2 + p t ■ p B ) 2 /E 2 - p 2 B , where p B is 



the B momentum and (Ei,pi) is the four-momentum of 
the initial state. The tties distribution for signal events 
peaks near the B mass with a resolution of 2.4MeV/c 2 , 
while the AE distribution peaks near zero with a reso- 
lution of 19McV. Events in the interval (-30, +60) MeV 
centered on —4.9 MeV in AE are accepted. The shift 
in the central value corresponds to the mean of the sig- 
nal AE distribution observed in data for the channel 
B + — > Z?°7r + , D° — > K + ir~ . An asymmetric interval is 
chosen to reduce the amount of BB background from 
4-body decays. We require events to lie in the range 
5.20 < toes < 5.29GeV/c 2 . This range is used for an ex- 
tended maximum likelihood fit to the toes distribution to 
determine the number of signal and background events in 
our data sample. The region is then subdivided into two 
areas: a sideband region (5.20 < toes < 5.26GeV/c 2 ) 
used to study the background Dalitz-plot distribution 
and the signal region (5.271 < to E s < 5.287 GeV/c 2 ) 
where the Dalitz-plot analysis is performed. Following 
the calculation of these kinematic variables the B can- 
didates are refitted with their mass constrained to the 
world average value of the £?-meson mass in order to 
improve the Dalitz-plot position resolution. 

The dominant source of background comes from light 
quark and charm continuum production (e + e~ — > qq). 
This background is suppressed by requirements on event- 
shape variables calculated in the T(AS) rest frame. For 
continuum background the distribution of |cos#t|, the 
cosine of the angle between the thrust axis of the se- 
lected B candidate and the thrust axis of the rest of the 
event, is strongly peaked towards unity whereas the dis- 
tribution is uniform for signal events. Additionally, we 
compute a Fisher discriminant T |19| , a linear combina- 
tion of five variables: Legendre polynomial moments Lq 
and L 2 , the absolute value of the cosine of the angle 
between the direction of the B and the detector axis mea- 
sured in the CM frame, the absolute value of the cosine 
of the angle between the B thrust axis and the detec- 
tor axis measured in the CM frame and the output of a 
multivariate i3-fiavor tagging algorithm [2l| . The Fisher 
coefficients are calculated from samples of off-resonance 
data and signal Monte Carlo (MC). The selection require- 
ments placed on |cos#t| and optimized with MC sim- 
ulated events, accept 51% of signal events while rejecting 
95% of background events. 

Other backgrounds arise from BB events. There are 
four main sources: combinatorial background from three 
unrelated tracks; three- and four-body B decays in- 
volving an intermediate D meson; charmless two- and 
four-body decays with an extra or missing particle and 
three-body decays with one or more particles misiden- 
tified. We veto candidates from charm and charmo- 
nium decays with large branching fractions by reject- 
ing events that have invariant masses in the ranges: 
2.97 < m 7r±7r = F < 3.17GeV/c 2 , 3.56 < m,i lT < 
3.76GeV/c 2 , 1.8 < m K ±^ < 1.9GeV/c 2 and 1.8 < 
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m,±,T < 1.9GeV/c 2 . These ranges reject decays 

from J/tp -> £+£~, %j){2S) -> and I} ^> X+tt- (or 

7r + 7r~) respectively, where € is a lepton that has been 
misidentified. 



We study the remaining charm backgrounds that es- 
cape the vetoes and the backgrounds from charmless B 
decays with a large sample of MC-simulated BB decays 
equivalent to approximately five times the integrated lu- 
minosity of the data sample. The 73 decay modes that 
have at least one event that passes the selection crite- 
ria are further studied with exclusive MC samples. Of 
these 54 are found to be significant, and the MC sam- 
ples of these modes are used to determine the Toes and 
Dalitz-plot distributions that are used in the likelihood 
fits. These distributions are normalized to the number 
of predicted events in the final data sample, which we 
estimate using the reconstruction efficiencies determined 
from the MC, the number of BB pairs in our data sample, 
and the branching fractions listed by the Particle Data 
Group 01 an d the Heavy Flavor Averaging Group [2^ . 
The predicted yields of BB background events in the 
signal region are 263 ± 16 (270 ± 16) for the negatively 
charged (positively charged) sample. 

To extract the signal and qq background fractions 
we perform an extended maximum likelihood fit to the 
toes distributions. The signal is modeled with a double 
Gaussian function. The parameters of this function are 
obtained from a sample of nonresonant K ± -k^-k ± MC 
events and are fixed except for the mean of the core 
Gaussian distribution, which is allowed to float. The 
qq toes distribution is modeled with the experimentally 
motivated ARGUS function The endpoint for this 
ARGUS function is fixed to \fsj2 but the parameter de- 
scribing the shape is left floating. The BB background 
toes distribution is modeled as the sum of an ARGUS 
function and a Gaussian distribution, whose parameters 
are obtained from the BB MC samples and are fixed in 
the fit. The yields of signal and qq events are allowed 
to float in the final fit to the data while the yield of BB 
background events is fixed to the value determined above. 

The results of the fits to both the negatively charged 
and positively charged samples are shown in Fig. ^ The 
fit yields 1047±56 (1078±56) signal events and 1016±25 
(999 ± 24) qq events for the negative (positive) sample in 
the signal region. 

We independently fit the Dalitz plot of the negatively 
charged and positively charged samples to extract the 
magnitude and phase of the intermediate resonances 
and the nonresonant contribution using an unbinned 
maximum-likelihood fit. We construct a likelihood func- 
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FIG. 1: The mes distribution, together with the fitted PDFs: 
the data are the black points with statistical error bars, the 
lower solid (red/dark) area is the qq component, the middle 
solid (green/light) area is the BB background contribution, 
while the upper blue line shows the total fit result. The up- 
per (lower) plot is for negatively charged (positively charged) 
events. 



tion: 



C(x,y) = (1 - f q g - f BB ) 



(13) 



EfciCje^Fj(a;,y) e(x,y) 



+ iW 



JV 



=1 U J C 



i Fj{x,y)\ 2 e(x,y) dxdy 



Q(x,y) 



J Q(x,y) dxdy 



+ !bb 



B(x,y) 



J B(x,y) dxdy'' 



where N is the number of resonant and nonresonant com- 
ponents in the model; e(x,y) is the signal reconstruction 
efficiency defined for all points in the Dalitz plot; Q(x, y) 
is the distribution of qq continuum background; B(x,y) 
is the distribution of BB background; and f q q and f BB 
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are the fractions of qq continuum and BB background 
components determined as described above and fixed in 
the Dalitz-plot fit. 

To allow comparison among experiments we present fit 
fractions (FF) rather than amplitude magnitudes where 
the fit fraction is defined as the integral of a single de- 
cay amplitude squared divided by the coherent matrix 
element squared for the complete Dalitz plot: 



30 r 



25- 



f 20-;; 

154; 

10 



J | Cje l6j Fj (x, y) | dxdy 



J V 'V* /•>•//; 



(14) 



5~ 



dxdy 



The sum of all the fit fractions is not necessarily unity due 
to the potential presence of net constructive or destruc- 
tive interference. The fit fraction asymmetry is defined 
as the difference over the sum of the B~ — > K~ir + ir~ and 
B + — > K + ir~ir + fit fractions: 
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3± 



FFj ~ FF, 
FF, + FF, 



(15) 



2: Background subtracted Dalitz plot of the combined 
B ± — > K if^ir data sample in the signal region. The plot 
shows bins with greater than zero entries, the area of the 
boxes being proportional to the number of entries. 



The qq continuum and BB backgrounds are modeled 



as two-dimensional histograms (m K ^ 



„- vs. mz 



the B + and B~ samples with bins of size 0.4( GeV/c 2 j" x 
0.4(GeV/c 2 ) 2 with linear interpolation applied between 
bins. The BB background distributions are taken from 
MC studies and the qq distribution is taken from the jties 
sideband data. We expect from MC studies 1515±59 BB 
background events in the sideband region (776 ± 37 for 
negative events and 739 ± 36 for positive events) , which 
is 10.8% of the reconstructed sideband events. The dis- 
tribution of these events is subtracted from the qq dis- 
tribution. In order to increase the statistical precision of 
this distribution we then add off- resonance data events. 
This increases the sample size by 1202 events (605 for 
negative events and 597 for positive events). The Dalitz M. 
plot of the data in the signal region after subtraction of 
the two background distributions can be seen in Fig. [21 

The two-dimensional efficiency distribution over the (m-0 — ^k-k) ~ zwoFo 

Dalitz plot e(x,y) is calculated with 1.3 million 

B ± — ► i\T ± 7r =F 7r =t: nonresonant MC events. All selection where cot 5b = ^h + h r Q- We have used the following val- 



For most resonant amplitudes the pole masses and 
widths are taken from the standard Particle Data Group 
tables 0. However, there are no consistent measure- 
ments for the coupling constants g^, qk and the pole 
mass m of the / (980) 00,111. We em P lo y a likcli " 
hood scanning technique in order to determine the best- 
fit values: g n = 0.11, gn = 0.36 and mo = 0.965 GeV/c 2 . 
The + component of the K + ir~ spectrum, which we 
denote (Ktt)q°, is poorly understood jH H<|; we 
use the LASS parameterization H3,|2^| which consists of 
the ^q°(1430) resonance together with an effective range 
nonresonant component. 



m K Tr 



q cot 5 1 



iq 



(16) 



+ e 



2iS E 



m r ^ 



criteria are applied except for those corresponding to 
the invariant mass veto regions. The quotient is taken 
of two histograms, the denominator containing the true 
Dalitz-plot distribution of the MC events and the numer- 
ator containing the reconstructed MC events with correc- 
tions applied for differences between MC and data in the 
particle identification and tracking efficiencies. The effi- 
ciency shows very little variation across the majority of 
the Dalitz plot but there are decreases towards the cor- 
ners where one of the particles has a low momentum. No 
difference in efficiency is seen between B + and B~ . The 
effect of experimental resolution on the signal model is 
neglected since the resonances under consideration are 
sufficiently broad. The average reconstruction efficiency 
for events in the signal box for the nonresonant MC sam- 
ple is 16.7%. 



ues for the scattering length and effective range param- 
eters of this distribution: a = 2.07 ± 0.10 (GeV/c) -1 and 
r = 3.32 ± 0.34 (GeV/c)" 1 It has been shown in the 
decay B — > J/ip Kir that the P — S phase behaviour well 
matches that observed in the LASS experiment [3(J. But 
since this parameterization is only tested up to around 
1.6 GeV/ c 2 we curtail the effective range term at the D° 
veto. Integrating separately the resonant part, the ef- 
fective range part and the coherent sum we find that 
the ifg°(1430) resonance accounts for 81%, the effective 
range term 45%, and the destructive interference between 
the two terms the excess 26% of (Ktt)q°. 

The nominal model comprises a phase-space non- 
resonant component and five intermediate resonance 
states: K * (892)tt + , (if7r)*°7r+, ^"(770)^+, f (980)K+, 
XcoK + . We choose this model using information from 
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previous studies [3l| and the change in the goodness-of-fit 
observed when omitting or adding resonances. The non- 
resonant component is modeled with a constant complex 
amplitude. Alternative models for the nonresonant com- 
ponents, such as that proposed in j^, were also tested 
and found to make negligible difference to the measured 
parameters. The results of the fit with the nominal six- 
component model are shown in Table[I]separately for B~ 
and B + data. The statistical errors are calculated from 
MC experiments where the events are generated from 
the PDFs used in the fit to data. The j o°(770) resonance 
shows the greatest difference in fit fractions between the 
two samples. The projection plots of the fit can be seen 
in Fig. 01 For the m^-±„.=F plots the requirement is made 
that m-TriTrT is greater than 2GeV/c 2 and vice versa in 
order to better illustrate the structures present. Using 
the fitted signal distribution we calculate the average re- 
construction efficiency for our signal sample to be 14.7%. 
This value, which includes all corrections due to differ- 
ences between data and MC, can be used to calculate the 
inclusive branching fraction from the signal yield results. 

As a measure of goodness-of-fit we evaluate the x 2 
across the Dalitz plot. A 75 by 75 bin histogram is used 
and a minimum of 10 entries per bin is required (for those 
cases where this requirement is not met then neighbour- 
ing bins are combined). The results for B~ (B + ) are a 
X 2 of 123 (148) for a total number of 116 (121) bins and 
10 free parameters in both cases. We also calculate the 
same x 2 neglecting the region 1.2 < m 7r ±^ T < 1.6 GeV/c 2 
where there is no resonant component in the nominal fit 
model. The results for B~ (B+) are a X 2 of 109 (112) 
for a total number of 108 (109) bins and again 10 free 
parameters in both cases. 

The omission of any of the nominal components from 
the fit results in a significantly worse negative log- 
likelihood with the fitted fractions and phases of the re- 
maining components varying outside their error bounds. 

We have tested for the sensitivity of these results 
to additional resonances that can be added in the 
fit function. In the 7T7T spectrum there are possible 
higher resonances including / 2 (1270), /o(1370), p°(1450), 
/o(1500)and /2(1525); in the Kir spectrum there are pos- 
sible i4T 2 *°(1430) and AT*°(1680) resonances. Each of these 
resonances is added in turn to the nominal signal model, 
to form an extended model, and the Dalitz-plot fit is re- 
peated. In general, adding another component does not 
significantly affect the measured fit fractions and phases 
of the six nominal components. We place upper limits on 
each of the possible additional components (Table [HJ. 

The systematic uncertainties that affect the measure- 
ment of the fit fractions and phases are evaluated sep- 
arately for B~ and B + . Each bin of the efficiency, qq 
and BB background histograms is fluctuated indepen- 
dently in accordance with its errors and the nominal fit 
repeated. The fractions of BB and qq events are var- 
ied in accordance with their errors and the fits repeated. 



TABLE I: Final results of fits, with statistical, systematic 
and model-dependence errors, to B~ and B + data with the 
six component model. 



Component 



Fit 



7T- 



Fit 



if*" (892) Fraction (%) 
/f"°(892) Magnitude 
A"* (892) Phase 
(A'jr)S" Fraction(%) 
(K-k)o° Magnitude 
(ifir)" Phase 



15.0 ± 1.6 ± 
1.0 FIXED 
0.0 FIXED 

49.fi ± 2.4 ± 0.7^0 



13.1 ±1.5 ±0.61 
1.0 FIXED 
0.0 FIXED 



0.01 ±0.12 ±0.03 ±0.06 



p"(770) Fraction (%) 
,3° (770) Magnitude 
p°(770) Phase 



/o(980) Fraction (%) 
/o(980) Magnitude 
/o(980) Phase 



10.5 ±1.7 ±0.41 



-0.55 ± 0.38 ± O.OSIS bs 
16.1 ±2.1 ±0.41'-; 



5.4 ±1.3 ±0.51 



1.19 ± 0.62 ± 0.161 



1.037 ± 0.080 ± 0.0261 
1.23 ±0.34 ±0.0712 



13.5 ± 1.9 ± 0.4i£J 
1.015 ± 0.095 ± O.OSOiiS 
2.29 ± 0.56 ±0.13in?« 



Xco Fraction (%) 
Xco Magnitude 
Xco Phase 



0.84 ±0.44 ±0.141 



2.55 ± 0.40 ± 0.09tjS 



1.20 ±0.50 ±0.131 



-2.52 ± 0.37 ± 0.3512;;? 



Nonresonant Fraction (%) 
Nonresonant Magnitude 
Nonresonant Phase 



4.3±1.3±0.8i;;j 
0.54 ±0.07 ±0.0612;} 



4.6 ±1.5 ±0.91^ 

0.60 ± o.ii ±0.06lS;S 



To confirm the fitting procedure, 500 MC experiments 
were performed in which the events arc generated from 
the PDFs used in the fit to data. A small fit bias is 
observed for the Xco and nonresonant components and 
is included in the systematic uncertainties. There is a 
contribution to the fit fraction asymmetry from possible 
detector charge bias, which has been estimated in previ- 
ous studies to be 1% |3Jj ■ A 0.4% systematic error on the 
total branching fraction comes from the error on the pre- 
dicted number of BB background events. The systematic 
uncertainties for particle identification and tracking effi- 
ciency corrections are 4.2% and 2.4% respectively. The 
calculation of N g g has a total uncertainty of 1.1%. The 
efficiency corrections due to the selection requirements on 
cos 8t, the Fisher discriminant, AE and tties have also 
been calculated from B + — > D Q ir + , D° — > K + ir~ data 
and MC samples, and the error on these corrections is 
incorporated into the branching fraction systematic un- 
certainties. 

In addition to the above systematic uncertainties we 
also calculate a mo del- dependence uncertainty that char- 
acterizes the uncertainty on the results due to elements of 
the signal Dalitz-plot model. The first of these elements 
consists of the two parameters of the LASS model of the 
Kit S-wave and is calculated by refitting adjusting the 
parameters of the LASS model within their experimental 
errors. The second consists of the three parameters of 
the /o(980) resonance and is evaluated by refitting with 
the para meter values measured by the BES collabora- 
tion [24| . The third element is due to the different possi- 
ble models for the nonresonant component and is evalu- 
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FIG. 3: Invariant mass projections for the data in the signal region and the fit results. The upper (lower) plots are for the 
B~ (B + ) sample. The left-hand (right-hand) plots show the m K ± lT (m^±^) spectrum. The data are the black points with 
statistical error bars, the lower solid (red/dark) histogram is the qq component, the middle solid (green/light) histogram is 
the BB background contribution, while the upper blue histogram shows the total fit result. The large dips in the spectra 
correspond to the charm vetoes. For the m K ± 7r ^ plots the requirement is made that m 7V ± 1T ^ is greater than 2GeV/c 2 and vice 
versa. 



ated by refitting; with the parameterization proposed by 
Belle |32j. The fourth element is the uncertainty due 
to the composition of the signal model and reflects ob- 
served changes in the parameters of the nominal compo- 
nents when the data are fitted with the extended models. 
The uncertainties from each of these elements are added 
in quadrature to give the final model-dependence uncer- 
tainty. 

In order to make comparisons with previous measure- 
ments and predictions from factorization models we mul- 
tiply each fit fraction by the total branching fraction to 
calculate the branching fraction of the mode. These 
branching fractions from each of the charge-separated 
fits are then averaged. For components that do not 
have statistically significant branching fractions 90% con- 
fidence level upper limits are determined. Upper lim- 
its are also calculated for the components added in the 
extended signal models. Upper limits are calculated 
from MC experiments where each experiment is gener- 



ated from the fitted PDFs but with all sources of sys- 
tematic uncertainty varied in accordance with their er- 
rors. The measured branching fractions, averaged over 
charge-conjugate states, and CP asymmetries from this 
analysis are summarized in Table [n] Charge conjugates 
are included implicitly throughout this table and the fol- 
lowing discussion. 

The total B ± -> K ± tt^tt branching fraction ((64.1 ± 
2.4 ±4.0) x 10 -6 ) has been measured with increased sta- 
tistical accuracy and is compatible with previous BABAR 
measurements. It differs from Belle's measurement of 
(46.6±2.1±4.3) x 10" 6 H3, which is significantly smaller, 
even after accounting for the XcoK mode, which Belle 
does not include. This result was cross checked by us- 
ing the same procedure to measure the B + — > D°ir + ; 
D° — > K + tt~ branching fraction, which was found to be 
consistent with the PDG value The total charge 

asymmetry has been measured to be consistent with zero 
to a higher degree of precision than previous measure- 
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TABLE II: Summary of measurements of branching fractions (averaged over charge conjugate states) and CP asymmetries. 
The first error is statistical, the second is systematic and the third represents the model dependence. 
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merits. 

After correcting for the secondary branch- 
ing fraction B(K*°(892) K+%-) = § we find 
the B+ -> if*°(892)7r+ branching fraction to be: 
(13.5 ± 1.2 ± 0.7i°J) x 10~ 6 . This is smaller than that 
measured in previous analyses that do not perform an 
amplitude fit to the Dalitz plot [3lll3^| but slightly larger 
than the value reported by Belle in their amplitude 
analysis [S2T |. The branching fraction measurement of 
B + — > p°(770)K + is the first measurement of the mode 
from BABAR and is consistent with that from Belle [32j |. 
It is also broadly compatible with many theoretical 
predictions H 0, S |I @. The 5+ ^ 70(980)^+ 
branchin g fr action is in good agreement with earlier 
analyses |3ll and the recent Belle amplitude analy- 
sis |32j . The forward-backward asymmetry apparent in 
both the K*°{892) and / (980) bands in Fig. is well 
reproduced by the fit and is not due to reconstruction 
efficiency effects but to S — P interference in the Dalitz 
plot. 

The (i^7r)o° component appears to be well modeled by 
the LASS parameterization, which consists of a nonreso- 
nant effective range term plus a relativistic Breit-Wigner 
term for the Kq°(1A30) resonance itself. Removing the 
phase-space nonresonant component from the nominal 
model gives very little change in the goodness-of-fit x 2 
or the fit likelihood. It is unclear whether this compo- 
nent is required in addition to the nonresonant part of 
the (Kiy)q° component. We can calculate the branching 
fraction for B + — > ifQ (1430)7r + using our knowledge of 
the composition of the (Kit)q component and find it 
to be: (44.4 ± 2.2 ± 2.0±^ ± 4.5) x 10~ 6 , where the 
fourth error is due to the uncertainty on the branching 
fraction of i4Tg°(1430) — * Kir combined with the uncer- 
tainty on the proportion of the (Ktt)q° component due to 
the Xq°(1430) resonance. The Belle collaboration finds a 



similarly large -fCo°(1430) branching fraction though they 
treat the Xq°(1430) as a separate component from the 
rest of the S-wave, modeled as a nonresonant component 
that has variation in magnitude but no variation in phase 
over the Dalitz plot j^- 

For B+ -> p°(77Q)K + and B+ -> X coK + the differ- 
ences in phase between the B + and B~ decays are 
1.74 ± 0.73 (2.4 standard deviations, a) and 1.21 ± 0.54 
(2.2<j) respectively, where the errors are statistical only. 
The if*°(892)7r + charge asymmetry is consistent with 
zero, as expected from the Standard Model predic- 
tions 0, E3 E3 IE El • There is no evidence of new physics 
entering the penguin diagram loop. 
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